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Abstract. The normality equations for the Newtonian dynamical systems on an 
arbitrary Riemannian manifold of the dimension n > 3 arc considered. Locally the 
solution of such equations reduces to three possible cases: in two of them the solution 
is written out explicitly, and in the third case the equations of normality are reduced 
to an ordinary differential equation of the second order. Some new examples of 
explicit solutions of normality equations are constructed. 



1. Introduction. 

The concept of dynamical systems that admit the normal shift was introduced in 
[1,2]. These are the systems which describe the motion of a mass point according to 
the Newton's second law mv ~ F(r,r) and which satisfy some specific geometrical 
condition that permits them to implement the normal shift of any hypersurface 
along their trajectories. Such geometrical condition was called the condition of 
normality. In [1,2] this geometrical condition of normality was brought to the 
form of a system of partial differential equations for the force field F(r,r) of the 
dynamical system. These equations were called the equations of normality. 

In [3,4] the equations of normality were generalized for the case of Newtonian 
dynamical systems on an arbitrary Riemannian manifold. Note that in this case 
the normality equations were written in covariant form. Therefore one can apply 
the methods of differential geometry to these equations. 

The geodesic flows of the metrics that is conformally equivalent to an original 
metric of the Riemannian manifold form a wide subclass of the dynamical systems 
accepting the normal shift. The problem of coincidence of trajectories of a dynam- 
ical system with the trajectories of such geodesic flow was called the problem of 
metrizahility. It was considered in [5,6] where an explicit form of the force field for 
all metrizable dynamical system that admit the normal shift was obtained. 

The symmetry analysis of the equation of normality was started in [7] where all 
classical symmetries of these equations in two-dimensional spatially homogeneous 
case were found and the corresponding self-similar (invariant) solutions were con- 
structed. As a result of the symmetry analysis some examples of non-trivial (non- 
metrizable) dynamical systems accepting the normal shift in two-dimensional case 
were found. The importance of the search for such dynamical systems especially for 
the dimension n > 3 was pointed out to authors by academician A. T. Fomcnko. 
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Some simplest examples of non-metrizable systems were; constructed in [8], however 
the systematical search for such systems has yet to be undertaken. 

The main goal of the present article is the complete local analysis of normality 
equations for the dimension n > 3 and the construction of new non-metrizable 
dynamical systems that admit the normal shift. 

2. The equations of normality. 

Let M be a n-dimensional Riemannian manifold with the metric tensor . Let's 
denote by .t^, . . . , the local coordinates of some point on M and by f ^, . . . , f" 
the coordinates of a tangent vector at this point. The equations of Newtonian 
dynamical system describing the motion of the point of unit mass m = 1 have the 
following form: 

(2.1) x'=v\ Vtv'=F\^,w). 

The normality equations for the force field F of the dynamical system (2.1) are 
written as two systems of differential equations: 

(2-2) I {ViFk+VkFi-2v-^FiFk)N^Pi+ 

, + v-^{VkFiF^ - WkF^N''NrFi)P^ = 



(2.3) 



P^VkF'iP^. = ^ . " P^ 



_ P^VkFiP ; 

« " n-1 

The systems of equations (2.2) and (2.3) were derived separately in [3] and [4]. 
The first system was called the equations of weak normality, and the second one 
was called the additional norm,ality equations. The equations of normality (2.2) and 

(2.3) are differential equations in the covariant derivatives V and V in the expanded 
algebra of tensor fields. A tensor field from the expanded algebra is distinguished 
from an ordinary tensor field on A4 by the doubled set of arguments: it depcnids on 
a point X € M as well as a tangent vector v at this point. Therefore there are two 
types of covariant derivatives in the expanded algebra: the spatial gradient V and 
the velocity gradient V. The calculation of components of the covariant derivative 
V for a tensor field U of (r, s)-type is simply the differentiation with respect to the 
components of velocity: 

. . . duy-Y 

(2.4) y^u;i:::;: = 

There are no connection components T^j in (2.4). Therefore the covariant deriva- 
tive V can be defined on a manifold without any metric. The components of the 
covariant derivative V are calculated in more complicated way: 
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The equations (2.2) and (2.3) arc the equations with respect to the force vector 
F. But apart from the components of the force vector, in (2.2) and (2.3) there are 
also the components of tensor fields N and P, and a scalar field v. These are given 
parameters in the normality equations: v is the field of velocity modulus v = |v|; 
N is the unit vector directed along v and P is the orthogonal projector on the 
hyperplane that is perpendicular to the vector v: 

(2.6) N' = v'/v, P] = 6} - N' Nj. 

Both covariant derivatives V and V of a metric tensor are equal to zero. This 
makes calculations in the expanded algebra similar to calculations in the ordinary 
algebra of tensor fields on M. The covariant derivatives of given fields v, N and P 
are calculated from the formulae (2.4) and (2.5). 

3. A SCALAR SUBSTITUTION. 

Consider the first equation from (2.2). Let us denote A = F'^Nj.. The scalar 
field A from the expanded algebra has the sense of the orthogonal projection of the 
force F on the direction of the velocity vector. So, the force field F can be written 
in the form: F = ^ N + F. Substituting this expansion into the first equation from 
(2.2), we find: 

(3.1) v-^FiPl + ViAPl=0. 

Here we've used the obvious equality P^ Ni — 0, which follows from (2.6), and the 
fact that N is a unit vector. From the orthogonality of the vectors F and N we have 
Fi Pg = Fq. So, taking into account everything told above, from (3.1) we obtain: 

(3.2) Fq = ANq-vPiViA. 

Prom (3.2) we see that all components of the force field F are defined by the scalar 
field A. 

Let's substitute (3.2) into the second equation of weak normality (2.2). It leads 
to quite cumbersome, but not complicated computations. As a result of these 
computations we obtain: 

Pi {ViA + V P''' VsAVk^iA- 
3.3 ^ , . . - - 

-vN''VkViA-N>'AVkViA) = 0. 

In order to derive (3.3) we've used the following relationships: 



(3.4) 




= 0, 




= Ni, 


(3.5) 




= 0, 




= Pt, 


(3.6) 




= 0, 


^iP^ 


= -v-^{Pig N^ + P^Nq) 



that are obtained by direct computations taking into account the relationships (2.4), 
(2.5) and (2.6). 
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Now let us substitute (3.2) into the additional equations of normality (2.3) and 
use once more the relationships (3.4), (3.5) and (3.6) in combination with Ni = 
and Pg — Pk- From the first equation (2.3) we derive 

_ AT'- VrVgA VkA - VkVqA) = 0. 
After the substitution (3.2) into the second equation (2.3) we obtain: 



P^^VkV,AP''^ = ^^^^^P^, 



(3.8) 

where the number n is a dimension of the manifold M. It's convenient to rewrite 
the last equation in the following form: 

(3.9) P^Vk^qAP'^^ =XPi, 

where A = A(x, v) is a function on the tangent bundle TM. It is the very form this 
equation was initially derived in [4]. 

4. Fiber spherical coordinates. 

The equations (3.3), (3.7) and (3.9) form the complete list of the equations of 
normality after the scalar substitution (3.2). The last equation from this list has an 
important difference from the other ones. If we write out the covariant derivatives 
VfcVgA, in explicit form, from (3.9) we obtain: 

'^''^™ = ^^«- 

The equations (4.1) contain only the derivatives of A with respect to the coordinates 
, . . . , i;" in the fiber of the tangent bundle TM . The fiber is a n-dimensional linear 
vector space with the metric gij which is constant within the fiber. 

Let's fix some point P on the manifold M . The conditions v = \v\ = const 
decompose the fiber of TM over P into the union of non-intersecting spheres S^~^ 
with the radii r = v. Let u^,. . . , m""^ be local coordinates on the unit sphere S^~^ 
and let v = N(u'^, . . . , m"~^) be the parametrical equation of the unit sphere in the 
fiber of TM over the point P. Then the relationships 

=u"- Ari(u\... ,u"-^), 

(4.2) 

= u".iV"(u\... ,tt"-i) 

define the change of Cartesian coordinates , . . . , v" for the spherical coordinates 
,. . . where w" = |v|. The choice of coordinates (4.2) in each fiber of TM 
defines the fiber spherical coordinates on the tangent bundle TM. The dependence 
of the choice of such coordinates upon the point P can be made smooth, however 
for analysis of the equations (4.1) this does not matter. Therefore we'll consider 
the relationships (4.2) taking the point P to be fixed. 
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Let's transform the equations (4.1) by making the change of variables (4.2). 
The metric Qij is constant within the fiber and the derivatives V, = d/dv^ are 
covariant derivatives in the metric Qij. The equations (4.1) are tensorial equations, 
therefore it's sufficient to recalculate the components of tensors P^, Pj, gij and the 
derivatives Vi in the new coordinates. 

Note that according to (4.2) we have dv'^/du" = Therefore vector N is the 
n-th coordinate vector in the spherical coordinates: N = (0, . . . ,0, 1). So, for the 
components of the projector P in spherical coordinates we have 

gin-l 

gn-ln-1 



in view ot wiiat we said above, m tlie splierical coordinates the equation (4.1) 
has the initial form (3.9). But the components of the connection in spherical 
coordinates w^, . . . are nonzero. For k in the range 1 < i,j,k < n — 1 the 
quantities do coincide with the components of the metric connection on the 
spheres |v| = const. The other components are calculated explicitly: 

By virtue of the structure of the matrices (4.3) for i = n or j = n the equations (3.9) 
are fulfilled identically. Taking into account (4.4) we can write the rest non-trivial 
equations as the equations in the covariant derivatives V on the spheres |v| = const: 

/I dA \ 

(4.5) v,V,.4=(^A+-— jffi,-. 

Since the function A is an indeterminate parameter in (4.5), it's convenient to 
introduce another parameter fj, and to write (4.5) more simply: 

(4.6) ^iVjA = figij. 

Consider a diff'erential consequence from the equation (4.6). Let us write (4.6) in 

the form VjVfeA = figjk then apply the operator Vi to both sides of this equation 
and alternate the obtained equation with respect to the indices i and j: 

(VjVj - VjVi)Vfe^ = Vingjk - Vjiigik- 

(4.7) - Rlij A = Vi/x gjk - ^jpi Qik- 

Here i?| = v~'^ {5f gjk — 5| gik) is the curvature tensor of the sphere S^~^ given 
by the equation ii = |v| = const. 

Let's contract the above equation (4.7) with the metric g^'^. It leads to the 
following relationship: 

n—2 - , , ^, - 
— 5- ViA = n - 2 Vi/x. 



(4.3) 



p. 





1 




911 

gn-l 1 





(4.4) 



^ nn 



0, 



an _ _9il 
-' V 
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Note that the additional equations of normality (2.3) arise in the dimension n > 3 
only, therefore in the obtained equations n — 2 ^ 0. Hence 

(4.8) -ViA = Vifi. 

Let's substitute the last expression for V,j4 into the equations (4.6). It leads to the 
following differential equations for the functional parameter fj,: 

(4.9) ViVj/z = -v~'^ Qij jJi. 

Theorem 4.1. The general solution of the system of equations (4.9) on the sphere 
|v| = const has the form: 

n 

(4.10) ti = ^miN\u\... ,u^), 

i=l 

where mi,... , m„ are some arbitrary constants and N^,... are the compo- 
nents of the vector N. 

Proof. Because the equations (4.9) and the expression (4.10) are covariant tensorial 
equalities it's convenient to check the statement of the theorem for some particular 
choice of coordinates u^, ... ,m" on the sphere |v| = const. Let's render concrete 
the change of variables (4.2), taking the metric gij to be brought to the unitary 
form gij = Sij in the coordinates v^,. . . , t;" in the fiber of TM over the considered 
fixed point P: 

= V sin((/,„,_i) siii(-u,„,_2) ■ • • • • sin(!i3) sin(M2) sin(Mi), 
v"^ = V sin(M„_i) sin(M„_2) • . . . • sin(u3) sin(u2) cos('Ui), 
= V sin(u„_i) sin(u„_2) • . . . • sin(w3) cos(m2), 
(4-11) 

ti"-i = V sin(M„_i) cos(u„_2), 

t;" = V COs(w„_i). 

In spherical coordinates u^, . . . , w" defined by the relationships (4.11) the metric 
gij is diagonal, moreover ga > and the quantities Hi = y/gu are called the Lame 
coefficients: 

Hi = V sin(u„_i) sin(M„_2) • . . . • sin(M3) sin(M2), 
H2 =v sin(M„_i) sin(u„_2) • . . . • sin(u3), 

(4.12) 

Hn-2 = V sm(u„_i), 

Hn-l = V, 
Hn =1. 

The components of the metrical connection are defined by the Lame coefficients 
(4.12) in the following way: 

Hi dHj ^ 

~ Hk du^ ^ Hk dui {Hkf du^ 
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Most of the components of the connection i?^^ are zero. The calculation of the other 
ones makes possible to write (4.9) explicitly. As a result we obtain the equations 
of two types. The equations of the first type correspond to the case i < j < n — 1 
in (4.9): 

/ . , „x cos(w,) 

^'■''^ 

Here /Xj and /Xjj are the partial derivatives of ^ of the first and the second order 
with respect to the corresponding variables and . The equations of the second 
type correspond to the case i = j < n — 1: 

(4-14) M»+ E 3^ n --^K) /^^=- n --^K) 

fe=i+l ^ ' \p=i+l / \p=i+l / 

Note that, when i = n — 1, the equation (4.14) is reduced to the form na = —fi. 

Let us make the further proof of the theorem by induction on n. Wc take, as the 
base of induction, the case n = 2. In the initial geometrical problem the additional 
equations of normality arise only for the dimension n > 2. Therefore the equations 
(4.6) are also absent. The relationship (4.8), which lets us bring (4.6) to (4.9), was 
also derived for n > 2. But once the equations (4.9) are already derived, their 
extrapolation for the case n = 2 has no contradiction. 

So, when n = 2, the systems (4.13) and (4.14) are reduced to the equation: 
/Uii = —ji. The general solution of this equation is: 

(4.15) = f^i sin('u^) + m2 cos(u^), 

where mi and m2 are constants. In the case n = 2 the relationships (4.11) have 
the form: 

(4.16) = V = V s.\n{ui) , v'^ = v = v cos{ui) . 

Comparison of (4.15) with (4.10) and taking into account (4.16) show that the 
statement of theorem 4.1 is fulfilled for the case n = 2. 

Now let n > 2. Consider the last equation from (4.14) with i = n—l. It has the 
form = —jj. Hence 

(4.17) M = A sin(u"~^) + m„ cos(u"~^), 

where jl and m„ are parameters that in general case depend on u^, . . . , w"~^. Let's 

substitute (4.17) into the equations (4.13), setting j = n — 1. As a result the con- 
tribution of the first summand from (4.17) is canceled and we obtain the following 
relationships for m„: 

„ " =0 for all i = 1, . . . ,n — 2. 

Hence m„ = const. Therefore the substitution of (4.17) into the rest equations 
(4.13) and (4.14) leads to the equations for the parameter jl, which have the same 
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form as (4.13) and (4.14). but with n — 1 in place of n. It means that one can apply 
the assumption of induction. As a result we obtain: 

(n-2 \ n-1 /ri-2 \ 

sm{uP) + ^Wfe sin(wf) cos{uP-^), 

p=l / fe=2 \p=k J 

where mi, . . . ,rn„_i are some constants. Now it remains to substitute (4.18) into 
(4.17) and to compare each of the summands in the obtained formula with (4.11). 
From this comparison we see that (4.18) and (4.10) are the same. So, the theorem 
is proved. □ 

Now let's return to the relationship (4.8) and write it in the form Vi(A + u^ /n) = 
0, from where we obtain: 

(4.19) A + fj, = const . 

Note that (4.19) and parameters rrii from (4.10) are constants within the separate 
spheres |v| = const only. Therefore, denoting bi = |v|, rui, we can formulate the 
following statement. 

Theorem 4.2. The scalar field ^(x, v) from the extended algebra of tensor fields 
satisfying the equations of normality (3.9) has the form: 

(4.20) A = a(x, |v|) + |v| 6i(x, |v|) N\ 

where the scalar field a and the covector field b from the expanded algebra depend 
on the velocity modulus v = \v\ in the fibers of the tangent bundle TM. 

5. The refinement of the scalar substitution. 

The obtained formula (4.20) completely defines the dependence of the function 
A upon velocity v. The dependence of A upon coordinates should be refined by the 
substitution of (4.20) into the equations (3.3) and (3.7). Introduce the following 
notations a' = da/dv and a" = d'^a/dv'^. It is not difficult to check that a' and 
a" are scalar fields from the expanded algebra of tensor fields and they depend on 
velocity modulus only. Moreover 

(5.1) Via = a' Ni, Via' = a" Ni. 

The derivatives = dbk/dv and 6^' = db'j^/dv define two covector fields b' and b" 
from the expanded algebra of tensor fields. These fields also depend on velocity 
modulus only and there are following relationships analogous to (5.1): 

(5.2) yibk = b'kNi, Wib'k = blN,. 

Substituting (4.20) into the equations of normality (3.3) and taking into account 
(5.1) and (5.2), we obtain: 



(5.3) 



P'g {Via + a'bi-b'ia + vN'- Vibr -vN"- Vrh) = 
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The equation (5.3) contains the covariant derivatives V^a and Vi^r- The formula 

(2.5) is written in the variables , . . . ,x",v^,. . . But for the fields a and b 
the natural variables are . . . ,x",v, where v = \v\ = ^ gkq v'' vi. While in 
the variables x^,. . . , x", u^, . . . , the velocity modulus v depends on a;^, . . . , a;" 
because gkq depends on , . . . , x" . Therefore 

„ da a' doka h „ o-' 

In view of concordance of the metric and connection two last summands in this 
expression are canceled. An analogous cancellation happens in the computation of 
Wibr- For this reason the covariant derivatives of the fields a and b are calculated 
as if they contain no dependence on v: 

(5.4) ViO = da/dx\ = dW/dx' - hu- 

By virtue of (5.4) the covariant derivatives V^a and Vj^r depend on the velocity 
modulus only. 

Let's replace the projector in (5.4) by the explicit expression = S^ — Nq 
from the formula (2.6) 

VqU + a'bq-b'ga + v {Vqbr - Vrbq) N'' 

^^'^^ - Nq{Vra + a'br-bla)N'' = 0. 

There are three groups of summands in the equation (5.5). The first three sum- 
mands depend on the velocity modulus only. The other ones depend on the modulus 
of the velocity vector and on its direction as well. The direction of v is given by 
N. One of these summands is linear in N, while another is quadratic with respect 
to N. The substitution of N by — N does not change the modulus of the velocity 
vector. This substitution changes a sign of the summand that is linear in N and 
the rest summands in (5.4) remain unaltered. Therefore the summand linear with 
respect to N vanishes identically. Hence 

(5.6) Vqbr - Vrbq = 0. 

The summand quadratic with respect to N contains the contraction of the vector 
N and the covector with the components V^a + a' br — b'^a which depend on |v| 
only. While the modulus of the velocity vector is unchanged, the vector N can be 
made orthogonal to this covector. So, the summand quadratic in N is also zero: 

(5.7) V,ia + a'bq-b'ga = 0. 

Conclusion: the equation of normality (3.3) for the field A of the form (4.20) is 
reduced to the equations (5.6) and (5.7) for the fields a and b. 

The substitution of (4.20) into (3.7) add nothing new. It leads to the equation 
that is the consequence of (5.6). For this reason (5.6) and (5.7) remain as unique 
basic equations of normality for the fields a and b. Taking into account the rela- 
tionships (5.4) and symmetry of the connection components T^j enables us to bring 
(5.6) to the form: 

The covector field b(x, v) from the expanded algebra of tensor fields depending only 
on the modulus of velocity can be treated as the one-parameter family of 1-forms 
on the manifold: h = bi dx^, where t; is a parameter. Then the equation (5.8) is a 
closedness condition for each form from this family. 
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Example 1. Let a = 0. Then the equation (5.7) is fulfilled for any choice of b. 
So, each one- parameter family of closed 1-forms corresponds to some Newtonian 
dynamical system accepting the normal shift. The force field of such system has 
the following form: 

(5.9) Fg = \v\bi{^,\v\)i2N'Nq-5i). 

Note that the simplest examples of the dynamical systems accepting the normal 
shift from the papers [1], [2] and [8] are systems of the form (5.9). 

Example 2. Let b = 0. Then the equation (5.7) has the form: V^a — 0. It means 
that the function a does not depend on a point of the manifold M, i.e. a = a{v). 
For the force field of the corresponding dynamical system we have: 

(5.10) F, = a'(|v|)7V,. 

This is a trivial class of the dynamical systems that admit the normal shift. The 
trajectories of the systems (5.10) coincide with geodesies of the manifold M. 

Now consider the case a ^ and b ^ 0. In this case the equation (5.7) is not 
trivial. Let's write it in the following form: 

(5-11) (^ + ^^l;)i^H = 6^. 

Let us set Xq = d/dx'^ + bq d/dv and compute the commutator of two such differ- 
ential operators, taking into accoimt the relationships (5.8): 

(5.12) [Xk, Xq] = {bk b'^ - bq b'k) d/dv. 

The equations (5.11) arc overdctcrmincd. Using (5.8) and (5.12), from (5.11) we 
derived the following differential consequence: 

(5.13) bjb''-^b'\=bq(bl-^b'^. 



Lemma 5.1. The components of two covectors b and c satisfy the relationships 
bk Cq = bq Ck if and only if they are linear dependent. 

Proof. If one of the covectors is equal to zero the lemma statement is trivial. Let 

b 7^ 0. Then, at least one component of the covector b is not zero. Suppose bi ^ 0. 
Then from bkCq = bqCk we derive Ck = Xbq, where A = ci/bi. It means c = Ab. 
So, the lemma is proved. □ 

Let us apply Lemma 5.1 to the equation (5.13) where b ^ 0. So, as a consequence 
of (5.13) we obtain the following equation: 

(5.14) b'^-^b'^=Xbq, 

where A is some scalar. Let's differentiate the equation (5.11) with respect to v and 
write the obtained result in the form: 



a^lnlol, a^ini 



(5-15) K-^K=^^^>>q + 



dip' dv dx'i 
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The comparison of (5.14) with (5.15) leads to tlic equation that is also to be con- 
sidered as a differential consequence of the equations (5.7) and (5.9): 

where the scalar fi is obtained from A by subtracting the second logarithmic deriv- 
ative of the function a with respect to v. 

For further analysis of the obtained equations (5.15) note that the equations 
(5.8) are locally solvable. Any field b satisfying the equations (5.8) is defined by 
some scalar field /3 from the expanded algebra which depends on the modulus of 
velocity only: 

(5.16) br = VrP = df3/dx''. 

The substitution (5.16) into the equation (5.15) transforms this equation to the 
form: 

/r-i^x 9 fa'\ dp 

Lemma 5.2. // spatial gradients of two functions a{x,v) and (3{x,v) are pro- 
portional: da/dx'^ = ^d(3/dx'^ then in some neighborhood of points at which the 
gradient of the function (3 is not zero, there is the representation a = F{(3, v) where 
F is some function of two variables. 

Proof. In Lemma 5.2 the quantity v plays the role of a parameter. Therefore to 
prove the lemma it is convenient to assume v by fixed w = wo = const. Then, 
while grad/3 ^ 0, the local coordinates .t^,... , on M can be taken so that 
x^ = (3{x,vo). In this case from the proportionahty of the gradients a and /? we 
have: 

da da ^ 

dx'^ dx"^ 

It means a = a{x^,VQ). And the hmctional dependence of a on a; -'^ and uq in 
these specially chosen coordinates defines the function F{P,v). The choice of such 
coordinates can be made smoothly depending on the parameter vq- Therefore 
F{l3, v) is the smooth function of two variables that defines a relation between a 
and /? in the form of following relationship: a — F((3,v). □ 

Let's apply Lemma 5.2 to the equations (5.17). As a result we obtain the differ- 
ential equation binding a and /?: 

(5.18) ^=^(/^'^)- 
Let's substitute (5.18) into the equation (5.11). Then 

^ln|a| ^.^ , dp d^p 

Let $(/?, v) be the antiderivative of the function F{P, v) with respect to P for fixed 
V. In other words, the function $ is connected with F by the relationship: 



(5.20) 
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that defines v) up to the summand depending on v only: 

(5.21) $(/?, v) $(/?, v) + *(w). 

After the substitution (5.20) into (5.19) taking into account that /3 = /3(x, u) lets 
us to rewrite the equation (5.19) in the following form: 

(5.22) ^L|«| + $(/?,,)_|^)=0. 



By virtue of (5.22) the expression in parentheses is the quantity depending on v 
only. Therefore because of the arbitrariness (5.21) in the choice of function $(/3, v) 
we have: 

(5.23) ln\a\ = ^-mv). 

The substitution of (5.23) into the equation (5.18) leads to the differential equation 
that for the given $(/3, v) defines the dependence of /3 on v: 

(5.24) ^-=55^(,- + l) + 2£gZ), 

The relation of the quantities f3 and $(/3,f) and the force field of the dynamical 
system contains an element of arbitrariness (see (5.21)). This arbitrariness defines 
the following gauge transformations that do not change the form of the equation 

(5.24) : 

/3^^(x,t;)=/3(x,i;)+V(^;), 

(5.25) _ . 

$ ^ $(/3, u) = $(/3 - V(v), u) + V/(f )• 
The transformations (5.25) do not change the quantities a = ln\a\ in (5.23). The 
equation (5.24) can be rewritten in the form of a system of two equations of the 
first order 

r p' = a + <i>{f3,v), 

(5.26) I , ^ dHP,v) 

I " dp ■ 

Note that the system (5.26) also admits the gauge transformations (5.25) comple- 
mented by the relationship a a{v) = a{v). 

In general case the equation (5.24) as well as the system (5.26) is not possible to 
solve explicitly. However this equation enables to characterize exactly the extent of 
arbitrariness in the definition of force field of the dynamical system accepting the 
normal shift. For fixed choice of the function $(/?, t;) the general solution of the 
equation (5.24) contains two integration constants / and h: 

(5.27) Piv)=B^{v,f,h). 

The parameters / and h in (5.27) can depend on the coordinates , . . . ,x". So, 
/ and h are two scalar fields on the manifold M: f — f[x) and h = /i(x). The 
corresponding force field of the dynamical system for (5.27) has the following com- 
ponents: 

The formula (5.28) is considerably less effective than the formulae (5.9) and (5.10) 
for the cases considered earlier. It contains the function i?$ which is not arbitrary. 
Some special cases when this formula can be made more effective we'll consider in 
following section. 
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Theorem 5.1. The force fi,eld F oj the dynamical system that admit the norm,al 
shift is locally defined by one of three formulae (5.9), (5.10) or (5.28). In the last 
case it contains three arbitrary parameters: two scalar fields /(x) and h{x) and the 
function of two parameters $(/3, v). 

Example 3. In [5] and [6] the class of mctrizablc dynamical systems that admit 
the normal shift was studied. Consider these systems in the context of the above 
construction. Let us choose the function $(/3, u) of the special form: 

(5.29) i>(/3,w) = -lni?(we-'^/"), 

where H = H{^) is some smooth function of one variable. The substitution of 

(5.29) into (5.24) leads to the equation: 

(5.30) ^" = ir'~''^ {^'-^)- 

Finding the general solution of the equation (5.30) in explicit form is too problem- 
atic. However one particular solution is easily guessed: 

(5.30) p{^,v)=vf{x). 

Here / is some scalar field on the manifold M. The substitution of (5.29) and (5.30) 
into (5.23) defines the field a. The field b is also easily calculated: 

(5.32) a = H{ve-^')ef , bq = vVJ. 

The fields (5.32) define dynamical systems with the force field of the form: 

F, = -|vr VJ + 2 Vkfv" Vi + Ng H{\v\ e"^') e^' . 

This is the force field of the mctrizable dynamical systems from [6] . It contains one 
arbitrary scalar field / and an arbitrary function H of one variable. It's clear that 
this case does not exhaust the functional arbitrariness declared by theorem 5.1. 
This fact confirms the existence of non-trivial non-metrizable dynamical systems 
accepting the normal shift in the dimension n > 3. 



6. Some new examples. 

Example 4. Consider the case, when the equation (5.24) becomes linear. So, let 
us choose the function ^{(3,v) being hncar with respect to /3 and set 

(6.1) $(/j,„) = _£M^ + ln^». 

The expression (6.1) is not a general form of the function linear in f3, however 
the general case can be reduced to (6.1) by the gauge transformation (5.25). The 
system of equations (5.6) equivalent to the equation (5.24) for the function ^{f3,v) 
of the form (6.1) is written as: 



(6.2) 
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The system (6.2) is integrated easily. As constants of integration two scalar fields 
/(x) and h{x) arise: 



(6.3) a = -ln0»+ln/., /?= ^"\y/^ . 

Hence for the force field F of the dynamical system by virtue of (4.20) and the 
scalar substitution (3.2) we obtain the expression: 

--^«.-m(^|^.^).».v,-.;). 

Consider the equation (5.24) once again. The general solution of this equation 
depends on two integration constants: /? = B^(v, f,h). Let the parameter h be 
functionally expressed via the parameter /, i.e. h = h{f). This diminishes the 
number of arbitrary parameters and leads to the function 

(6.5) = 0{vJ) = B^{vJ,h{f)), 

that defines some one-parameter subfamily in two-parameter family of solutions 
of the equation (5.24). Let's consider (6.5) as a function of two variables and 
differentiate it with respect to v. As a result we obtain one more function of two 
variables: /?' = (3'{v,f). Without loss of generality, one can suppose that the 
dependence on / in (6.5) is locally reversible: / = f{(3,v). Let's substitute this 
into the function /?' = (3'{v,f). As a result we obtain the following differential 
equation of the first order: 



(6.6) l3' = U{l3,v), 

where U{P,v) = fi' {v, f{(5,v)). The equation (6.6) is compatible with (5.24). One- 
parameter family of its solutions is exactly the subfamily (6.5) of solutions of the 
equation (5.24). Let us differentiate (6.6) with respect to v and substitute the result 
into (5.24). It leads to the relationship: 

(6.7) «^ UiP, V) + ^ = ^ m0, .) + !)+ ^^(^' ^) 



df3 ' dv 8(3 ^ ^"^^ ' ' dv ' 

that binds the functions U{(3,v) and $(/3, t;). The relationship (6.7) is the compat- 
ibility condition of two ordinary differential equations (6.6) and (5.24). 

Note that (6.7) is one equation for two functions. To solve (6.7) let's introduce a 
new function W{(3,v) = U{l3,v) — $(/?,«). Then the equation (6.7) can be written 
in the form: 

Let's make the following substitution into the differential equation (6.8): 



(6.9) W{l3,v) = -lnw{l3,v), 
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After this substitution the equation (6.8) becomes quite simple: 

du dw dw 

The general solution of the equation (6.10) is defined locally by one arbitrary func- 
tion of two variables (ji(/3, i)): 



dl3 dv' dp' 

Hence for the function U{P,v) and for the function W{P,v) introduced above we 
obtain: 

Let us write the first relationship from (6.11) in the form of an equation for the 
function ^(/3, v): 

Such equation is solved by the methods of characteristics. The characteristics of 
the equation (6.12) are defined by the following ordinary differential equation 

(6.13) /3' = UiP,v) + l, 

and the function (p{p,v) by virtue of (6.12) is the first integral (conservation law) 

of the equation (6.13). Hence a value of the hmction cj) can be used to parameterize 
the family of solutions of the equation (6.13). Let's write it as: 

(6.14) P = B{v,(p). 

The function (6.14) differs from (6.5). It inverts the functional dependence </> on /3, 
i.e. B{v,(f){p,v)) = p. Differentiating this relationship we obtain: 

(6 15) dBdl^^ ^ + ^ ^ ^ 

dcf) 8(3 ' dv d(j) dv 

Let us express the derivatives dcp/dp and dcp/dv from (6.15) via dB/dcj) and dB/dv 
and substitute the result into (6.11). It gives: 

Form (6.16) it is easy to define the function $ but in the variables v and (j): 
(6.17) ^ -W = B^ -\n{B4,) -I. 

In the same variables v and (j) it is naturally to write the equation (6.6) too. So, we 
obtain the equation defining the dependence of ^ on t;: 

B^[v, 0) 

It is derived from (j)' = 8(j)/d(3-U + d4>/dv and from the relationships (6.15) defining 
partial derivatives (p with respect to v and 0. 

Summarizing the computations just done, note that the initial equation (5.24) 
contains an arbitrary function $(/?, ti), we've expressed it via (j){P,v) and then 
via the function B{v,(j)). Hence it is convenient to consider B{v,(j)) as a primary 
function and to express everything in terms of this function. 



16 



ANDREY YU. BOLDIN, RUSLAN A. SHARIPOV 



Example 5. Let us render concrete the choice of the function B{v, (p) to solve the 
equation (6.18) explicitly. Let 

B{v,(l)) = In 



3v \(t) + 3 

Then the equation (6.18) for the function (j){v) gets the form of the Bernoulli equa- 
tion: 

The general solution of this equation contains one integration constant / that de- 
pends on the space variables x^,. . . , and defines a scalar field on the manifold 
M: 

'^ = '^^"'^) = exp(-3^V2-2/)-l- 

Now it remains to compute the fields a and b in (4.20) to define the field A of the 
scalar substitution (3.2) 

a = e- = S, = -i-sinh^ (-^ - / 

' dx' d(j) df 3v 

Now the expression for the components of force field of the dynamical system ac- 
cepting the normal shift is easily derived: 



4 . ^2f 3|v|2 Vif{2N'Ng-5i) 

smh 

9 V 



(6.19) Fg = -— sinh^ ( ^ - / ) iV, 



The force field (6.19) contains one scalar field / only and this fact distinguishes this 

force field from the general case (5.28) and from the case (6.4). It happens because 
we've replaced the equation of the second order (5.24) by the equation of the first 
order (6.6). 



8. Final remarks. 

Three cases defined by the formulae (5.9), (5.10) and (5.28) characterize com- 
pletely the local construction of the dynamical systems accepting the normal shift 
on Riemannian manifolds of the dimension n > 3. Two dimensional case is dis- 
tinguished by a larger extent of diversity. In multidimensional case the system of 
normality equations is strongly overdetermined (however it is compatible). Hence 
the set of its solutions is more restricted. 

Note that the existence and number of local solutions of normality equations are 
not connected with the features of metric: the constantness of curvature, the num- 
ber of isometrics etc. Hence the existence and number of smooth global solutions 
of these equations depend exclusively on the topology of the manifold M. 
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